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ABSTRACT 

We show that there is a very simple relationship between differential and dimensional 
renormalization of low-order Feynman graphs in renormalizable massless quantum field the- 
ories. The beauty of the differential approach is that it achieves the same finite results as 
dimensional renormalization without the need to modify the space time dimension. 



The method of "differential regularization and renormalization" was introduced in Ref. [1] 
by Freedman, Johnson and Latorre (FJL) as a novel technique for regulating the ultraviolet 
divergences of quantum field theories and computing finite renormalized correlation functions. 
The method is applied to real space (rather than momentum space) Feynman diagrams, 
and requires no explicit cut-off and no explicit counterterms. Renormalized amplitudes are 
calculated directly, and these satisfy Callan-Symanzik-type equations which determine the 
renormalization group functions (3(g) and j(g) [1,2]. More recent investigations have explored 
the compatibility with supersymmetry [3] and gauge invariance [4], the extension to massive 
quantum field theories [5] , and the issue of lower-dimensional gauge theories [6] . 

In this Letter we explain how differential regularization and renormalization is related 
to dimensional regularization and renormalization. Such an understanding is clearly impor- 
tant since dimensional regularization and renormalization is the best understood and most 
widely-used procedure for analyzing renormalizable quantum field theories (particularly those 
involving gauge symmetry). Furthermore, a major motivation [1,2] for the development of 
differential renormalization is the observation that (at least formally) the method is compat- 
ible with both gauge and chiral symmetry, and may therefore prove more convenient than 
dimensional renormalization for chiral gauge theories as well as dimension-specific theories 
such as Chern-Simons, Wess-Zumino-Witten, etc. 

We shall illustrate the relationship between differential and dimensional regularization 
and renormalization in three types of renormalizable massless quantum field theories: 

(i) massless scalar theory (c/) 4 in four dimensions, </> 3 in six dimensions, and 4> G in three 
dimensions) ; 

(ii) pure Yang-Mills theory in four-dimensions, using the background field method; and 
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(iii) massless four-dimensional quantum electrodynamics. 

We shall restrict our attention to low-order graphs in order to pinpoint most clearly the 
essential correspondence between differential and dimensional renormalization (furthermore, 
much of the initial "processing" of higher-order graphs is common to both approaches, and 
has been discussed in detail in Refs. [1,2]). We shall also clarify the issue of gauge invariance 
in massless QED, where differential renormalization introduces (in principle separate) mass 
scales for the vertex and self-energy graphs which must however be correlated in a specific 
way to satisfy the one-loop Ward identity. This is to be contrasted with the dimensional 
renormalization approach in which there is a single mass scale \i (arising from the fact that 
the coupling constant acquires a non-zero engineering dimension) and the Ward identity is 
satisfied automatically (in, for example, the minimal subtraction scheme). 

The essential idea [1] of differential regularization and renormalization stems from the 
observation that most primitively divergent low-order Feynman graphs are well-defined in real 
space except for short-distance singularities at coincident points. If these singularities are too 
severe, the graph will not have a well-defined Fourier transform and higher-order graphs of 
which the original graph is a sub-graph will contain divergent integrals — hence the need for 
some form of regularization. The strategy of differential regularization and renormalization 
is to isolate the coincident point singularities and to reduce the degree of divergence of these 
singularities by expressing the singular terms as derivatives of less singular terms. In practice, 
this strategy is implemented in two distinct stages. 

In the first stage no dimensionful regulator is needed. In this initial stage one may use 
straightforward algebraic manipulations and simple differential identities in order to isolate 
the coincident point singularities and to extract derivatives so as to reduce the degree of 



divergence. In massless quantum field theories (in d 7^ 2) we typically encounter power-law- 
type singularities in bare graphs, and these may have their degree of divergence reduced by 
using the identity 

□I l-P+2 
T 
1 r /I \ 

W = 7 ^V73 V » C 1 ) 



(-p + 2)(d-p) 

In general (and, as we shall see, in particular for renormalizable theories) one confront the 
singularity |a;| _d , whose degree of divergence may only be reduced further by the introduction 
of an arbitrary (but essential) logarithmic mass scale. *^ It is clear that one cannot use Eq. (1) 
as it stands because of the ^j— pole. 

In the second stage of differential regularization and renormalization one proceeds by 
defining the regulated form: 
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1 (\iyM 2 \x\ 2 
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This relation is certainly true away from the origin, and, furthermore, the dependence on the 
arbitrary (but necessary) mass scale M yields the appropriate 5-function singularity at the 
origin: 
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To obtain Eq. (3) we have used the fact that the massless scalar Green's function in d (7^ 2) 
dimensions is 



* Note also that \x\~ d is the borderline between powers of |x| which have or do not have a 
well-defined Fourier transform (see below, Eq. (7)). 

' In a super-renormalizable theory this may not be necessary — it may be possible to apply 
differential renormalization without introducing a mass scale. See Ref. [6] for an example. 



and it satisfies 

UG{x) =5 {d) {x) . (5) 

It is worth stressing that in differential regularization and renormalization there are no infinite 
(or finite) counterterms — one may simply replace \x\~ d by its regulated form as in (2). Then 
the scaling relation (3) is the key to deriving the Callan-Symanzik equations satisfied by the 
renormalization amplitudes [1] . In the differential approach, one also integrates by parts freely 
within graphs and in taking Fourier transforms [1]. The consistency of this procedure and 
its compatibility with unitarity has been explicitly checked (in massless </> 4 theory through 
three-loop order) in a cutoff version of differential renormalization [7] . 

This two-stage implementation of the differential renormalization strategy has been ap- 
plied with impressive success and efficiency to many examples of renormalizable quantum 
field theories [1-7]. We now show that this same strategy may be implemented using real- 
space dimensional regularization and renormalization. Recall that real-space dimensional 
renormalization of quantum field theories, known as the "method of uniqueness," [8,9,10] 
is an enormously powerful calculational tool with which multi-loop computations may be 
performed with far greater ease than with conventional (principally momentum space) tech- 
niques. Here, however, our emphasis and motivation are quite different --we are interested 
in the comparison with differential renormalization, and consequently we are interested in not 
only the infinite parts of graphs, but also their regulated finite parts. 

The first important point to realize is that the initial stage of "processing" a bare Feyn- 
man graph (as outlined above) does not require a regulating mass scale, and may be performed 
in exactly the same manner in an arbitrary dimension d of space-time. It is only in the second 
stage, when one encounters the singularity |x| _d , that any difference between dimensional and 



differential regularization arises. The dimensional regularization approach described here is 
simply another way of regulating this singularity. In fact, if the theory is formulated in the 
non-integer space-time dimension d = D — 2e (where D is the "real" integer dimension and 
e is infinitesimal) then the initial processing of the bare graph does not lead to a singularity 
\x\~ d but to a singularity |x| _d+re , where re is some integer multiple of e. Indeed, since in 
a massless d = (D — 2e) -dimensional theory the coupling constant acquires a mass dimen- 
sion proportional to e, one in fact encounters a singularity fj, re \x\~ +re where // is a universal 
dimensionful parameter associated with the coupling constant.* Now one is free to use the 
identity (1) directly, isolating the ^— pole as a 1/e pole: 

1 1 

// re M _d + re — —n re ni-T-i _d + re + 2 

e r(2 — d + re) 



r (2 - d + re) T 



(d_ 1 ) [ ' 2(2-d) LJ { \x\ d -* J + 



Thus, in the dimensional approach, the x = singularity of the left-hand side of Eq. (6) 
is regulated as a "counterterm" 5 < -- d \x) : with both a 1/e pole and a finite e° coefficient, 
and a finite regulated e° term of exactly the same form as the differential regularization 
expression (2) (with the dimensional regularization mass scale (i identified with the differential 
regularization mass scale M). As a matter of terminology, we shall refer to the 1/e pole 
counterterm as the "infinite counterterm," the e° counterterm as the "finite counterterm," 
and the e° regulated term involving \i as the "non-counterterm finite part." This is precisely 
analogous to the situation in conventional momentum space dimensional regularization, where 
the divergences of divergent momentum integrals are regulated as Laurent expansions in e. 



* The fact that the "correct" power of \i appears with the bare graph is a special feature 
of renormalizable theories. 
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Indeed, we may equivalently present Eq. (6) in momentum space form by Fourier transforming 
the left-hand side: 



2-7r d / 2 r (- N i / x >< 



d d xe ik - x /j r£ \x\- d+re = - \2J ( V 



, d — re\ \ \k\ 
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(7) 



Here -0(2) = -^ lnT(z) is the digamma function and 7^; = —ip(l) = 0.5772 ... is Euler's 
constant. This agrees with the Fourier transform of the right-hand side of Eq. (6) computed 
by freely integrating by parts (as is usual in dimensional regularization) . In fact, using Eq. (7) 
it is straightforward to convert any of the real space expressions we discuss in this paper into 
the corresponding momentum space formulas. 

Renormalization now proceeds exactly as usual in the dimensional regularization ap- 
proach (albeit in real space rather than in momentum space) - - i.e. by defining an appropri- 
ate subtraction scheme (see e.g. Ref. [11]). Note from Eq. (3) that the scaling dependence of 
the "non-counterterm finite part" produces contributions of the "finite counterterm" form. In 
more conventional language this simply corresponds to adjusting the renormalization scale. 

There are,of course, minor technical complications (see below) in gauge theories involv- 
ing some tensor algebra and Lie algebra manipulations and in fermionic theories involving 
some Dirac matrix algebra, but the essence of the correspondence between differential and 
dimensional regularization (at low order) lies in Eqs. (2), (3) and (6). 
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Massless Scalar Theories 



Consider the Euclidean space massless scalar theory with action 



s = / d''.v { !<£□< 



n\ 



(8) 



where d 7^ 2 and n^2. The scalar Green's function is as in (4) and (5) and so, apart from 
tadpole contributions, the lowest order two-point function is (see Fig. 1) 
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For relevant values of d and n, this graph is singular at the origin and it may have its degree 
of divergence reduced by applying Eq. (1). For the theory to be renormalizable (at this order) 
this must lead to a \38(x) counterterm — this requirement relates the space-time dimension 
d and the interaction power n as: 



(2-d)(2-n) =4 



(10) 



(where we have used that D|x| 2_d ~ S^ d '(x)). Thus, for example, we have (/> 4 theory in 
four dimensions, <p 3 in six dimensions or (p 6 in three dimensions. In these cases, the two- 
point function is proportional to |x| _2_d . Applying formula (1) once, we arrive at a \x\~ d 
singularity. Then using (2) we may simply write down the finite renormalized amplitude in 
the differential approach: 
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(11) 



In the dimensional approach, with d = D — 2e (and the integer part, D, of the space-time 
dimension satisfies (10)) the coupling constant A is replaced by A — *■ Ao/U*™ -2 -* 6 , where Ao is 
dimensionless and \x is a universal mass scale.* Then 
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Using the relation (6) and expanding in powers of e we obtain 
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Thus we see that the dimensional renormalization two-point function agrees with the 
differential renormalization result if we choose our dimensional renormalization scheme to be 
one of subtracting both the infinite and finite counterterms, leaving just the non-counterterm 
finite part, with the scales M and \x identified. And because of the scaling relation (3), any 
other scheme which retains any part of the finite counterterm simply corresponds to rescaling 
the arbitrary mass scale M. 

2. Four-dimensional Yang— Mills Theory in Background Field Method 

In Ref. [1], Freedman, Johnson and Latorre (FJL) used differential regularization and 
renormalization to compute the one-loop two-point function in the background field method. 



* Another way to regulate this massless \<p n theory is to keep d at its integer value but 
take n to be non-integer, n = N + 6 with N an integer: this is the "5-expansion" of Bender 
and collaborators [12]. 
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Recall [13], that due to the explicit background field gauge invariance of the background 
field method, the two-point function is sufficient to determine the renormalization group (3- 
function. The one-loop contribution to the effective action from the gluon and ghost loops in 
Fig. 2 is (after performing the Lie algebra contractions) 



Q 2 (B) = - [(2.a) + (2.b)] 



x 



9 ^-^C A I d*xd d yB a Jx)B«(i,) 



d 



(14) 



--d„d u (G(x - y) 2 ) + d (d^Gix - y)) (d u G(x - y)) 
+ Ad^d v (G(x - y) 2 ) - 45 F n (G(x - yf) 

Here B^ are the background gauge field potentials,* Ca is the quadratic Casimir in the 
adjoint representation, and all derivatives are with respect to x. G(x — y) is the massless 
scalar Green's function as in Eq. (4). After performing some simple tensor algebra in (14) we 
find the manifestly transverse expression 
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Now, setting d = 4 — 2e and then using the relation (6) once we obtain 
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Note that B a coincides with A a of Ref. [13], while in Ref. [1] B a = gA a 
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The 1/e pole part of (16) (from which one may deduce the one-loop /3-function) agrees with the 
momentum space dimensional regularization computation [13], and the non-counterterm finite 
part (also from which one may deduce the one-loop /3-function) agrees with the differential 
regularization result of FJL (see Eq. (II.G.22) of Ref. [1]), with M = \x. Once again, in terms 
of dimensional renormalization, differential renormalization consists of a subtraction scheme 
in which only the non-counterterm finite part is retained as the renormalized amplitude. 

3. Massless Four-Dimensional QED 

In d dimensions, the massless fermion Green's function is 

SW = ~ (^) G{x) (17) 

where G(x) is the massless scalar Green's function in (4), and so S(x) satisfies 

1 .JLs(x) = UG(x)= S^ ( x ) . (18) 

ox 

Then the one-loop self-energy graph in d = 4 — 2e dimensions is (see Fig. 3) 

E(x) = - (e // e ) 7pS[x)'y v 6pvG(x) . (19) 

Using the Dirac matrix relation 7^71,7^ — (d — 2)7^, we find 



2. .2en ,\ f r ( 1 ~ € )\ _, d |_|-4+4e 



S(x)=e ^ £ (l-e)( ^ 2 _ e \ 1-q^x 



en /I A d e , . e 2 . <9 _ /lnn 2 |:d 2 \ „, . 

i^U + ^ +1+ln ^ 7 '^ (x) -^ 7 ^ n (^p ) +G(e) 



(20) 



16tt 
where in the last step we have used the relation (6). This agrees with the standard 

{e.g. Ref. [11]) one-loop momentum space dimensional regularization result, and the non- 
counterterm finite piece agrees with the differential renormalization result of FJL [1], with 
the differential renormalization mass scale M^ taken equal to \i. 
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The analysis of the one-loop QED vertex function (see Fig. 4) in real space requires 
some preliminary algebraic and tensor manipulations in order to isolate the short-distance 
singularities. (The reason for this is simply that the vertex effectively depends on two space- 
time coordinates, rather than just one as in all the previous examples.) However, these "stage 
one" manipulations are exactly the same as in the differential regularization approach (see 
[1]), except that we work in d rather than four dimensions. The vertex is (see Fig. 4) 

o 

Vj(x, y, z) = (e y 2 ~ d/2 j 7 M 5'(x - z)-y 3 S(z - yhJ^G(x - y) 

(ft \ / ft \ \ ) 

-q^G(x - z)\ (g^G(z - y)\ G(x - y) , 

where in the last line we have specialized to d = 4 — 2e dimensions and we have used the 
Dirac matrix identity: ^^aljlbl^ = ^Ibljla - (4 - rf)7 a 7j76- 
Using the notation of FJL, we write the vertex as 

Vj(x, y, z) = -el/j 3e (2767j-7a - 2e7 a 7 i 7 6 ) V ab (x -z,y-z) , (22) 

thereby defining 



Vab(u, v) ee ( — G(«) ) ( — b G(v) ) G{u - v) 



Exactly as in FJL, one then isolates the singularity in V ab within the trace part by writing 
(for example) 



A (Gin) ( JL G („)) G(u - „)) - ^ ( G („) G („) (± 



( d d S ab d d\ 

S ab d d 

+ -G(u)G( V )-.-G(u- V ) 

V ab (u,v)+ S -^G(u)G(v)^.-^G(u-v) . 
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Note [1] that V ab (u,v) has a finite Fourier transform in four dimensions, so we only need 
to regulate the remainder, V ab — V ab . Using the defining relation (5) for the scalar Green's 
function we find 



\ab 



v ab - v ab 



(4-26) V 4tt 2 



r(l-e) 



S(u — v)\u\ 



-4+4e 



(24) 



In the differential approach one would have e = and so one would proceed (see [1]) by using 
the regulated form of |zx| 4 given by Eq. (2). In the dimensional regularization approach we 
expand \u\~ + as in Eq. (6) (with a factor jU 2e "borrowed" from the full vertex expression in 
(22); remembering that the vertex has an overall mass factor n e due to the mass dimension 
of the coupling e = eo/U e ) to obtain 
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If we set d = 4, then V* b (u, v) | is the differential regulated form of V ab found by FJL (see 



I reg 



Eq. (II. C. 9) of Ref. [1]). Inserting (25) into the full vertex expression (22) we find 



V J (x,y,z) = V* JL (x,y,z) 



+ e 0/ u 



reg 



e e 
167T 
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JE + - +ln7T 



lj 8{x-z)8{z-y) , (26) 



where K FJL I is the final differential regulated form in [11, with the differential renormal- 

J I reg ° L J ' 

ization mass scale for the vertex, My, taken equal to \x. (Notice that the [7^ + | + ln7r] e° 
coefficient in (25) becomes [7^ + \ + ln7r] in (26) because of the contraction of the - pole in 
V ab with the e^aljlb prefactor in the full vertex expression in (22).) Once again we see that 
the differential regularization and renormalization approach yields the non-counterterm finite 
part of the dimensional regularization approach. 
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It is straightforward to check that the dimensional renormalization vertex in (26) and 
self-energy in (20) satisfy the one-loop Ward identity 

d i V j (x,y,z)=(6(z-x)-6(z-y))Z(x-y) . (27) 



dzJ 



This is immediately clear for the - pole part of the counterterms, but takes a little work for 
the other terms. Indeed, at first sight, comparing the [7^ + \ + In7r] factor in (26) with the 
"corresponding" [7^ + 1 + hi7r] factor in (20), it looks as though there is a discrepancy in 
the finite part. However, differentiating K? JL produces an extra finite counterterm con- 

* ' ° J I reg 

tribution which "corrects" matters. In contrast, in the differential renormalization approach 
there are no counterterms at all in E or Vj, so this extra finite "counterterm" contribution 
to -fp-Vj must be cancelled. This may be achieved simply by rescaling the masses - - if 
In (My,/ My) = \ the one-loop Ward identity is satisfied [1]. 

The implication of this is simply that while dimensional regularization and renormal- 
ization (for example, with a pole subtraction renormalization prescription) automatically 
respects gauge invariance (i.e. the Ward identity is satisfied), the differential regularization 
and renormalization procedure (which has no counterterms at all) will, in general, require 
correlations between the mass scales appearing in various regulated graphs in order to pre- 
serve gauge invariance. This is not in the least surprising, as it is clear from Eqs. (2) and (3) 
that adjusting the differential renormalization mass scale in a renormalizable theory produces 
finite counterterms. 

To conclude, we have shown that there is a very simple relationship between differential 
and dimensional renormalization in low order graphs. It is not unreasonable to expect that a 
systematic application of these ideas to higher order graphs may yield a proof that differential 
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renormalization works just as well as dimensional renormalization to all orders. However, the 
more interesting task is to apply differential renormalization to theories in which dimensional 
regularization is problematical. 
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FIGURE CAPTIONS 

Fig. 1: Lowest order two-point function in a \(p n massless scalar field theory. There are (n — 1) 
internal lines. 

Fig. 2: One- loop graphs involving two external background fields. Wavy lines indicate quantum 
gluons and dashed lines indicate ghosts. 

Fig. 3: One-loop self-energy graph in QED. Straight lines indicate fermion propagators and wavy 
lines indicate photon propagators. 

Fig. 4: One-loop vertex graph in QED. The index j is the space-time index of the vertex. 
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